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An Overview of Integral Quadratic Constraints for 
Delayed Nonlinear and Parameter-Varying Systems 

Harald Pfifer and Peter Seiler 


Abstract —A general framework is presented for analyzing the 
stability and performance of nonlinear and linear parameter 
varying (LPV) time delayed systems. First, the input/output 
behavior of the time delay operator is bonnded in the frequency 
domain by integral quadratic constraints (IQCs). A constant 
delay is a linear, time-invariant system and this leads to a simple, 
intuitive interpretation for these frequency domain constraints. 
This simple interpretation is nsed to derive new IQCs for 
both constant and varying delays. Second, the performance of 
nonlinear and LPV delayed systems is bounded using dissipation 
inequalities that Incorporate IQCs. This step makes use of recent 
results that show, under mild technical conditions, that an IQC 
has an equivalent representation as a Unite-horizon time-domain 
constraint. Numerical examples are provided to demonstrate the 
effectiveness of the method for both class of systems. 


I. Introduction 

This paper presents a framework to analyze nonlinear or 
linear parameter varying (LPV) time-delayed systems. In this 
framework the system is separated into a nonlinear or LPV 
system in feedback with a time delay. Stability and perfor¬ 
mance is considered for both constant and varying delays. 
The analysis uses the concept of integral quadratic constraints 
(IQCs) ll26l . Specifically, IQCs describe the behavior of a 
system in the frequency domain in terms of an integral 
constraint on the Fourier transforms of the input/output signals. 
Several IQCs valid for constant and varying delays have 
already appeared in the literature, see e.g. m, ESI, ESI. 

This paper has two main contributions. The first contribution 
is to provide a simple interpretation for IQCs used to describe 
constant time delays. In particular, constant time delays are 
linear, time-invariant (LTI) systems and hence they have an 
equivalent frequency response representation. Thus IQCs valid 
for constant delays can, in most cases, be interpreted as a 
frequency dependent circle in the Nyquist plane. It is noted that 
this interpretation previously appeared in the robust control 
literature: “G”-scales for robustness analysis with real param¬ 
eter uncertainty can be interpreted with circles in the complex 
plane Q, 13. Here, the geometric interpretation is used to 
construct a new IQC valid for constant delays. Moreover, the 
frequency domain intepretation provides insight for generating 
a new IQC valid for time-varying delays even though such 
varying delays are not LTI. All these results are contained in 
Section |lV] 

The second contribution of this paper is to apply general 
IQCs for analysis of nonlinear and LPV delayed systems. 
The standard IQC stability theorem in 12^ was formulated 
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with frequency domain conditions and hence requires the 
“nominal” part of the feedback interconnection to be an 
LTI system. An application of this stability theorem to LTI 
systems in feedback with a constant delay is given in El- 
Previous work on delayed nonlinear systems bounded the 
nonlinear elements of the system and the time delays by IQCs 
and considered this frequency domain approach to analyze 
a “nominal” LTI systems under IQCs, see e.g. llMll . Here, 
dissipation inequality conditions are derived to assess the 
stability and performance of “nominal” nonlinear and LPV 
systems in feedback with a delay. The dissipation inequalities 
are time-domain conditions but IQCs are typically expressed as 
frequency domain constraints. Thus the key technical issue is 
that the analysis approach requires an equivalent time-domain 
interpretation for an IQC. Previous work along these lines for 
constant IQCs has appeared in Chapter 8 of M- In fact, a 
large class of IQCs, under mild technical conditions, have an 
equivalent expression as a finite-horizon, time-domain integral 
as recently shown in Ell, ms. This time-domain expression 
enables IQCs to be easily incorporated into a dissipation 
inequality condition as shown in Section [V] These analysis 
conditions can be formulated and efficiently solved as sum-of- 
squares optimizations and semidefinite programs ||4| for 
nonlinear and LPV delayed systems, respectively. Numerical 
examples for both system types are given in Section VI These 
results complement recent robust performance conditions for 
LPV systems ll38l . 11211 . ||3^ . This paper builds upon ll35l . In 
addition to the results in 051 . it includes a detailed description 
of the generation of IQC multipliers for time delays. 

There is a large body of literature on time-delayed systems 
as summarized in ESI, Q. The most closely related work 
is that contained in E3l, ESI, El, E3 which use IQCs 
to derive stability conditions for LTI systems with constant 
or varying delays. As noted above, the contribution of this 
paper is to extend these results to nonlinear and LPV de¬ 
layed systems. Lyapunov theory is an alternative framework 
in the literature of time-delayed systems ESI, El, ES- 
This approach essentially constructs Lyapunov-Krasovskii or 
Lyapunov-Razumikhin functionals to assess the convergence 
of the free (initial-condition) response of the delayed sys¬ 
tem. Stability conditions for nonlinear E3, ED and LPV 
i3 delayed systems have been developed in the Lyapunov 
framework. These methods treat the time delay as integrated 
with the dynamics of the plant. This is in contrast to the 
approach considered here which uses input-output stability 
(forced response) and separates the time delay from the 
“nominal” plant dynamics. A benefit of the IQC framework 
is that extends naturally to systems with many delays and/or 
uncertainties. On the other hand, some stability conditions 
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in the Lyapunov-Rrasovskii framework appear to use time- 
varying quadratic constraints that do not have counterparts in 
the IQC literature. This connection is not pursued here but 
may lead to new insights within the IQC framework. 


II. Notation 


K and C denote the set of real and complex numbers, re¬ 
spectively. KLoo denotes the set of rational functions with real 
coefficients that are proper and have no poles on the imaginary 
axis. KHoo is the subset of functions in KLoo that are analytic 
in the closed right half of the complex plane. C™^", 
and denote the sets of mxn matrices whose elements 

are in K, C, KLoo, KHoo, respectively. A single superscript 
index is used for vectors, e.g. M" denotes the set of n x 1 
vectors whose elements are in K. For z G C, z denotes 
the complex conjugate of z. For a matrix M G C™^", 
denotes the transpose and Af* denotes the complex conjugate 
transpose. The para-Hermitian conjugate of G € 
denoted as G'^, is defined by G'^(s) := G(—s)*. Note that on 
the imaginary axis, G'^(jtv) = G{juj)*. L2 '^he space 

of functions v : [0, 00) —>■ K" satisfying ||u|| < 00 where 


v{t)'^v{t) dt 


LJo 


0.5 


( 1 ) 


Given v G L 2 [ 0 ,oo), vt denotes the truncated function: 


VT{t) := 


v{t) for t < r 
0 for t > r 


( 2 ) 


The extended space, denoted L2e, is the set of functions v such 
that Vt G L2 for all T > 0 . Finally, the Fourier Transform 
V := fF{v) maps the time domain signal v G L^lOjOo) to the 
frequency domain by 

poo 

v{juj) := / e~^^^v{t)dt ( 3 ) 

Jo 


III. Problem Formulation 

Consider the time-delay system given by the feedback 
interconnection of a nonlinear, time-invariant system G and a 
(constant) delay Vr as shown in Fig. The delay w = 

is defined by w{t) = v{t — r) where r specifies the delay. 
The input/output signals can, in general, be vector-valued with 
w{t), v{t) G K"’'. The remaining signals in the interconnection 
have dimensions d{t) G and e{t) G M"'. The feedback 
interconnection is obtained by closing the upper channels of 
G with the time delay Vt- This feedback interconnection, 
denoted as Fu{G,'Dt), gives a time-delay system with input 
d and output e. 

A robust stability approach to analysis is pursued in this 
paper. Thus it will be more convenient to express the system in 
terms of the deviation between the delayed and the (nominal) 
undelayed signal, 5t(u) := T’t(i') — v. A loop transforma¬ 
tion, shown in Fig. can be used to express the feedback 
interconnection as F’„(G,iSr). This loop-shift amounts to the 
replacement w = w + v where w := Sr{v). The system G 


d 



e 


Fig. 1. Feedback interconnection with time delay 


obtained after this loop-shift is assumed to be described by 
the following finite-dimensional differential equation: 

XG = f{xG,W,d) 

v = hi{xG,w,d) ( 4 ) 

e = h2{xG,w,d) 

where xcit) G is the state of G at time t. 


d 



e 


Fig. 2. Loop transformation to F„{G,St) 

An input-output approach is used in this paper to analyze the 
stability and performance of the time-delay system. For a given 
delay r, the induced L2 gain for the feedback interconnection 
from d to e is defined as: 

||F,(G, 5^)11 := sup M (5) 

0^d&L^'‘[0,oo), a:G(0)=0 

It is important to note that the restriction to time t > 0 
implicitly assumes zero initial conditions for both Vt and St- 
Specifically, w = Vt{v) is more precisely defined on L2[0, 00) 
by w(t) = 0 for t G [ 0 ,r) and w(t) = v{t — r) for t > t. 
Similarly, w = St{v) is defined on ^2(0,00) by w(t) = —v{t) 
for t G [ 0 , r) and w(t) = v(t — t) — v(t) for t > t. The notion 
of finite gain stability used in this paper is defined next. 

Definition 1 . The feedback interconnection of G and St is 
stable if the interconnection is well-posed and if the mapping 
from d to e has finite L2 gain, i.e. there exists a finite constant 
7 > 0 such that ||F„(G, 5 r)|| < 7. 

Two main analysis problems are considered. First, determine 
the largest value of r such that the feedback interconnection 
is stable for all t G [ 0 ,t]. The first problem gives the delay 
margin for the system. Second, given a delay r less than the 
delay margin, determine the largest induced L2 gain from d 
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to e. The second problem gives the performance of the system 
for a fixed level of delay as measured by the L2 gain. 

Since the delay is constant, Vr defines a linear, time- 
invariant (LTI) system. In this case the delay has a well- 
known frequency domain representation. For constant delays 
w = can be expressed in the frequency domain as 

w{juj) = 'Dr{juj)v{joj) where 'Dr{joj) := In other 

words, the delay is equivalent to a frequency-by-frequency 
multiplication by 2 ?^. Similarly, has the frequency response 
Sr{juj) = — 1 . These frequency domain relations are 

used in the next section to derive simple, geometric constraints 
satisfied by the input/output signals of Sr- Additional technical 
details on the frequency domain can be found in standard 
textbooks, e.g. i). 

Up to this point the presentation has focused entirely 
on constant delays. However, the analysis problems can be 
extended to consider time-varying delays. The time-varying 
delay w = T>f,r{v) is defined by w{t) = v{t—T{t)) where T{t) 
specifies the delay at time t. The subscripts f and r denote that 
the delay satisfies r(f) S [0,f] and |T(f)| < r for all f > 0. 
In other words, f is the maximum delay and r bounds the 
rate of variation. If r = 0 then Vf^r corresponds to a constant 
delay with value r € [ 0 , f]. In addition, define w = Sf,riv) by 
w = T>r^r(v) — V, i.e. Sf^r is the deviation from the nominal 
undelayed signal. A time-varying delay is not a time-invariant 
system. Hence it does not have a valid frequency-response 
interpretation. However, the frequency-domain intuition can 
be used to derive constraints on the input/output signals of a 
time-varying delay (Section |IV-C|l. 


IV. Frequency Domain Inequalities 

This section describes different frequency domain con¬ 
straints on the delay operator that can be incorporated into 
the input/output analysis. The basic idea builds on common 
frequency domain inequalities that have appeared in the liter¬ 
ature, e.g. ED, ca, M- Given a constant delay r, define 
the LTI system (p as: 


Pis) := 2 


(sr)^ -I- 3.5sr -f 10 ® 
(st)2 -p 4.5st -1-7.1 


( 6 ) 


Fig. 3 (a) shows the Bode magnitude plots for Sr (solid 
line) and p (dashed line) The weight p is chosen to satisfy 
|5T0'a;)| < \(j>{juj)\ for all ui and hence Sr is a member of 
the following frequency weighted uncertainty set: 


{A : |A(ja;)| < \4i{juj)\ Vw} 


(7) 


This magnitude bound has simple geometric and algebraic 
interpretations at each frequency. The geometric interpretation 
of this uncertainty s et is a circle in the complex Nyquist plane 
as depicted in Fig. 3 (b) Sr{joj) follows a circle centered at 


— 1 with radius 1 (dashed circle). At each frequency Sr{joj) 
lies within the shaded circle of radius \(p{juj)\ centered at 
the origin. An algebraic interpretation can be given in terms 
of a quadratic constraint. The Fourier transforms for any 


input/output pair w = Sr{v) must satisfy the following 
quadratic inequality at each frequency: 


v{juj) 

* 


0 


1- 

1_ 

1 

_1 


0 

-1 


w[juj)_ 


This quadratic constraint is just a restatement of the norm 
bound on Sr at each frequency. Using this basic intuition, 
additional geometric constraints in the Nyquist plane can be 
expressed in the form of quadratic constraints at each fre¬ 
quency. Pointwise quadratic constraints are discussed further 
in Section |IV-A| for general LTI systems and the application 
to constant time delays is given in Section |IV-B| 



(a) Bode Magnitude Plot for Sr and bound (p 



(b) Circle Interpretation for |iST(fu;)| < |0(ya;)| 


Fig. 3. Norm Bound on Sr 


A. Pointwise Quadratic Constraints 

This section describes pointwise quadratic constraints for 
an LTI system A. For simplicity assume A is a single 
input, single output (SISO) system. The input/output rela¬ 
tion w = Av is represented in the frequency domain by 
w{juj) = A{juj)v{juj). This representation can be used to 
bound the input/output signals using frequency-by-frequency 
quadratic constraints (QC). 
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Definition 2 . Let 11 : jM —)■ be a Hermitian-valued 

function, called a ’’multiplier”. Two signals v,w G L2[0, c») 
satisfy the QC defined by If if the following inequality holds 


v{juj) 

w{juj) 


nO'w) 


v{jLj) 

w{juj) 


> 0, Vw 


(9) 


Moreover, the LTI system A satisfies the QC defined by If, 
denoted A G QC(n), if (|^ holds for all v G L2[0, oo) and 
w = Av. 


These QCs have an intuitive geometric interpretation. 
Specifically, the QCs can be interpreted as circle or half¬ 
plane constraints on the Nyquist plot of A{juj). Partition the 
2 x 2 multiplier as If = [ j. The diagonal entries are 

real since If is Hermitian. In addition, the QC is unaffected 
by positive scaling. Specifically, let A : jM. —M be a 
frequency dependent scaling that satisfies X{juj) > 0 Vui. Then 
A G QC(n) if and only if A G QC(An). Thus, without 
loss of generality, the multiplier can be normalized to have 
7 >'22(jtu) = — 1 , 0 , or +1 for all uj. This normalization provides 
a clearer geometric interpretation for the QC as described in 
the next lemma. 

Lemma 1 . Let A be an LTI system satisfying A G QC(n). 
At each frequency the QC can be normalized to one of three 
cases: 

1 ) If 7r22(ju;) = —1 then 

|A(ja;) - 7 r 2 i(ja;)| < x/tthQw) + | 7 r 2 r(jw)P ( 10 ) 

2 ) If 7r22(Juj) = 0 then 

0 < TTnO'uj) + A(juj)* 7 r 2 i(juj) + 7 r 2 i(ju>)A(Ju’) 

( 11 ) 


parameter uncertainty can also be interpreted as a circle 
constraint Q, a. 

Multiple QCs can be combined to obtain new QCs. If the 
LTI system A satisfies the QCs defined by then 

A also satisfies the QC defined by n(A) := X]/c=r 
for any real, non-negative numbers {Afc}^]^. n(A) is called 
a conic combination of the multipliers { 11 ^}^^. This fact 
enables many QCs on A to be incorporated into an analysis. 
However, it is important to recognize that such conic combi¬ 
nations have certain limitations. The main limitation is that 
any conic combination is, by Lemma just a circle or half¬ 
plane constraint in the complex plane. As a concrete example, 
consider the circle constraint on Sr described by Equation 
and shown in Fig. [ 3 ( 


The corTesponding multiplier 


IS 


Hi (jo;) := . Sr also satisfies the QC dehned by 

n2 = l}]. This second multiplier corresponds to a circle 

centered at —1 with unit radius, shown as the dashed circle 
in Fig. | 3 (b)| Thus Sr lies in the intersection of the shaded 
circle centered at the origin (defined by Hi) and the dashed 
circle centered at —1 (defined by 112). However, the conic 
combinations n(A) := AiHi -f A2n2 correspond to circles 
that “cover”, i.e. outer bound, this intersection. 


Remark. In fact, conic combinations of two QC multipliers 
are “tight”. Roughly, the conic combinations AiHi + A2n2 
define the smallest regions that contain the intersection of the 
sets defined by Hi and Il2. This statement is a geometric 
consequence of the S-procedure lossless theorem for complex 
constraints cni, a, uni. However the S-procedure is not 
lossless, in general, for three or more QCs. Hence conic 
combinations need not provide a tight bound on the set 
described by the intersection of three or more multipliers. 


3 ) If 7r22(iw) = -fl then 

|A(ja;) - 7 r 2 i(ja;)| > a/tthQw) -f | 7 r 2 i(jw)P (12) 


Proof. The QC in Equation must hold for all input signals 
V G L2[0, oo). Since w{juj) = A{juj)v{ju}), the QC can be 
rewritten as (dropping the dependence on joj): 


0 < 


1 

* 

1 


H 


A 


A 


— TTii -f A*7r2i + 7r2]^A — A*7r22A (13) 


If = 0 then this simplifies to Equation 11 If 


'^22{j<^) = ±1 then complete the square to express Equa¬ 
tion [T^ as in Equation [T^ or Equation □ 


Equation dehnes a circle and its interior (a disk) in 
the complex plane centered at 7r2i ( iuj) with radius given by 
•^/ttiiQ'o;) -I- |7r2i(ja;)p. Equation lljJ defines a half plane in 
the complex plane. For example if 7r2i(ja;) = 1 then 0 
defines the half plane given by Re(A{ja;)) > —^tthQu;). 
Finally, the case where 7122(^0;) = -|-1 corresponds to the 
non-convex set described by a circle and its exterior. Thus a 
QC defines a circle or half-plane constraint in the complex 
(Nyquist) plane at each frequency. The dissipation theory 
developed below is only valid for 7722 (jw) < 0 and hence 
the focus will be on quadratic constraints that are described 
by disks. “G”-scales used for robustness analysis with real 


The QCs defined above for SISO systems can be extended, 
with only notational changes, to multiple-input, multiple out¬ 
put (MIMO) systems. It will be sufficiently general for the 
time-delay analysis to consider repeated systems. Let A be 
a SISO, LTI system and dehne the n x n repeated system 
w = {A-In){v) by Wi = Avi for i = 1 ,... ,n. If A G QC(n) 
for a 2 X 2 multiplier H then the following QC holds for all 
V G £2(0,00) and w = (A ■ In)v: 


1 

* 

7711 (jw) • In 

7712 (jw) • In 

1- 

1_ 

w{juj)_ 


_772l(jw) ■ In 

7722 (jw) • In_ 



(14) 


Moreover A is LTI and hence A •/„ commutes with any n x 
n, frequency-dependent matrix D, i.e. D{A-In) = (A •/„)£). 
Thus the frequency-scaled system A := DAD~^ also satisfies 
the QC in Equation 14 Let (v,w) be any input-output pair 
for the scaled system A as shown in Fig. The associated 
input/output pair for the original system w = (A ■ In)v is 
related to the input/output pair for the scaled system hy w = 
Dw and v = Dv. Hence, A • /„ also satisfies the QC with any 
multiplier n„ : jR -G of the form 


n„(jw) 


7 rii(jw) • X{juj) TTi 2 {jU}) ■ X{juj) 
7721 (j W) ■ X{juj) 7722 (juj) ■ X (juj) 


(15) 
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where X{juj) := D{juj)*D{joj) > 0 . This more general, 
frequency-scaled multiplier can be used to reduce the con¬ 
servatism in the analysis. The use of X is analogous to 
the multipliers used in classical robustness analysis, e.g. the 
structured singular value p, Ell, 0, ED, na. 



A 



Fig. 4. Scaling of the Operator A ■ 


Fig. 5. “Small” Circle Constraint on St described by QCilla) 


B. Application to Constant Time Delays 

There are numerous circle constraints that can be used to 
bound St- A few examples are provided below to illustrate 
how the frequency-domain geometric interpretation can be 
used to impose constraints on St- The QCs on Sr can be 
converted, if needed, into equivalent QCs on Dr by reversing 
the loop-transformation, i.e. by replacing w = w — v m the 
QC. For clarity, the QC multipliers are given assuming Sr is 
SISO. However, as described previously, frequency-dependent 
(matrix) scalings can be introduced if Sr is MIMO. 

Example 1 . The Nyquist plot for St follows a circle centered 
at —1 with radius 1 . Hence, the most basic QC to bound St 
describes exactly this circle. This corresponds to the following 
multiplier; 

Hi "^1 . (16) 

-1 -1 


Example 2 . The multiplier in Example 1 does not depend 
on the value of the time delay t. Thus the multiplier in 
describes a very conservative constraint due to this delay 
independence. A simple delay dependent constraint is obtained 
from the frequency response 5 r(ja;) = — 1 . Thus at each 

frequency ST{juj) lies within a circle centered at the origin of 
radius | 5 r(ja;)|. The multiplier n2 for this circle is given by; 


112 (jw) 



2 

0 

0 

-1 


( 17 ) 


Example 3 . A smaller circle constraint can be constructed 
for St- The midpoint of the segment connecting ST{juj) and 
the origin is given by ^StUuj)- The following multiplier n3 
defines a circle centered at this midpoint with radius equal to 
the absolute value of this midpoint. 


naO’w) 


0 ^StUuj) 

^StUuj) -1 


The QC described by Hs is shown in Fig. 


( 18 ) 


As can be seen by the examples, there exist numerous ways 
to bound St using a QC. More examples are given in literature. 


For instance, in flCj a multiplier is provided that corresponds 
to a circle whose center moves along the imaginary axis 
with frequency. Unfortunately there are no general rules for 
which specific QC will provide the most useful stability and 
performance analysis results. It might seem intuitive to use the 
QC that describes the smallest area circle. However, this does 
not, in general, provide the least conservative analysis results. 
The best solution is to specify many different QCs and allow 
a numerical algorithm select the best conic combination. This 
will be described further in the next section. 

Two practical issues must be addressed to make these 
QCs useful for numerical analysis. First, recall that a system 
Fu{G,St) has delay margin of f if it is stable for all delays 
T G [ 0 ,f]. Thus the delay margin analysis requires a QC 
that covers St for all r G [ 0 ,f]. For example, the following 
multiplier covers St for all constant delays t G [ 0 ,f]; 


HaO'w) 


✓ 


2 -I 

0 


0 

-1 




( 19 ) 


4 0 

0 -1 


else 


At each low frequency (w < ^), II2 describes a circle 
centered at the origin that covers all Nyquist curves 
for r G [ 0 ,f]. This is essentially equivalent to the multiplier 
in Example 2 at low frequencies. However, these circles fail 
to have the desired covering property at higher frequencies. 
For w > 7 , the multiplier 112 is set equal to a circle centered 
at the origin of radius 2 . This ensures that the QC defined by 
n2 covers the entire Nyquist curve of St at high frequencies. 
The other QC multipliers specified in the examples can be 
similarly modified at high frequencies for use in delay margin 
analysis. 

The second practical issue is the need to approximate a QC 
multiplier with a rational function so that state-space numerical 
methods can be applied. For example 112 is a non-rational 
multiplier that describes a circle at each frequency. Fet (j)2{s) 
be any stable FTI system that satisfies |^20w)| A 
for w < 7 and \(p2{juj)\ > 2 for w > 7 . The specific choice in 
Equation|^satisfies these constraints. Then S^ satisfies the QC 
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defined by the following (rational) multiplier for all t G [ 0 , r]; 


112 (jw) 


l</>2(jw)P 0 

0 -1 


( 20 ) 


112 (jw) describes a circle with larger radius than 112 (jw) due 
to the choice of </> 2 (s). Hence Sr G QC(n 2 ) C QC(n 2 ) for 
all T G [0,f]. n 2 is the weighted multiplier (Equation]^ 
discussed earlier in the section. It is commonly used in the 
literature to analyze time delayed systems im, i26i, Qa. 


Remark. Numerical tools can be used to aid the construction 
of rational function approximations for more complicated 
multipliers. Briefly consider a (non-rational) multiplier H = 
^-1 ■ Lemmaj^the QC associated with this multiplier 
dehnes a circle centered at 7 r 2 i(ja;) with radius given by 
y/nii^ju}) + | 7 r 2 i(ja;)p. A rational approximation fl can be 
computed by the following procedure. Fit 7 r 2 i with a rational 
approximation 7 f 2 i, e-g- using the Matlab function fitfrd 
or some other numerical optimization. This approximation will 
introduce some fitting error and hence the radius of the rational 
multiplier may need to be increased to ensure Sr is covered. 
This would require tth to satisfy a lower bound constraint on 
its magnitude. A rational approximation with a lower bound 
magnitude constraint can be computed for tth, e.g. using 
the Matlab function fitmagfrd. Using this procedure, the 
following rational ht was constructed for Ha that covers Sr 
for all delays t G [0, f] 


naO'w) := 


0 

(jisijuj) -1 


where 


, ^ -2.19 (^) +9.02 (^)+0.089 

■= , - \ 2 

(^) - 5.64 (^) - 17.0 


( 21 ) 


( 22 ) 


Note that 7 r 2 i does not need to be stable, as is the case with 
the given choice of ^ 3 . 


C. Application to Varying Time Delays 

The QCs given for constant time delays hold at each 
frequency. These point-wise QCs fall within a more general 
framework based on integral quadratic constraints (IQCs) 
introduced in ES\ . A definition is now given for IQCs that 
extends the one given in Definition for QCs. 


Definition 3. Let U : jR ^ ^ 

Hermitian-valued function. Two signals v G L^^[0,oo) and 
w G [0, 00 ) satisfy the integral quadratic constraint (IQC) 
defined by H if 



v{juj) 

w(juj) 


n(juj) 


v(juj) 

w{jui) 


duj >0 


(23) 


where v{juj) and w{juj) are Fourier transforms of v and w, 
respectively. A bounded, causal operator A : 00 ) —>■ 

[0, c») satisfies the IQC defined by H, denoted A G 
IQC(n), if holds for all v G Lp [0, 00 ) and w = A(v). 


Clearly, if the QC holds pointwise in frequency it also 
holds when integrated over all frequencies, i.e. A G QC(n) 
implies A G IQC(n). The converse is not true and the 
more general IQC theory can be applied to nonlinear and/or 
time varying perturbations. IQCs were introduced in ll26l to 
provide a general framework for robustness analysis. Here, the 
focus will be on the use of IQCs to describe the input/output 
behavior of time-varying delays. A time-varying delay is not 
a time-invariant system. Hence the point-wise QCs cannot 
be used and the more general IQCs are required. However, 
the frequency domain arguments used to construct QCs for 
constant delays can provide intuition regarding IQCs for time- 
varying delays. The remainder of this section reviews known 
IQCs for time-varying delays ifTSll . llT9l . ifT^ . In addition a 
new IQC for time-varying delays is constructed using the 
frequency domain intuition gained from constant-delays. 

Recall the notation Df^r and Sr^r introduced for varying 
delays where f and r are bounds on the maximum delay and 
rate of variation, respectively. The basic IQCs for time-varying 
delays arise from two simple norm bounds. First, if r < 1 then 
the induced L 2 gain of the varying-time delay can be bounded 
by ||D^,r-|| ^ (Section 3.2 in ifTbll and Lemma 1 in ifT^ ). 
Second, let Sf,r o - denote Sr,r composed with an integrator 
at the input. The induced L 2 gain of this combined system 
can be bounded by ||5 t .rO 2 || < r (Lemma 1 in lITSl l. These 
two bounds are tight in the sense that the gain is achieved for 
some input signal v and time-varying delay r(f) that satisfies 
the bounds f and r (Lemma 1 in lfT9l ). Three IQCs are now 
provided for time-varying delays. These examples parallel the 
QCs provided in the previous subsection for constant delays. 
For clarity the multipliers are expressed assuming Dr^r and 
Sf.r are SISO. The extension to the MIMO case is discussed 
below. 


Example 4. Let w = Df^v) be a time-varying delay 
satisfying r < 1. The bound implies 

that llwll < ^^2— ||z;|| for all input signals v G L 2 [ 0 ,oo). 
After performing the loop-transformation w := Sf,r(,v), i.e. 
w = w + V, this inequality can be written as: 


poo 

v{t) 

T 

p- -1 

1 — r 


1 - 

■40 

1 _ 

Jo 



-1 -1 




(24) 


By Parseval’s theorem, this inequality can be equivalently 
expressed in the frequency domain as: 


poo 

v{jio) 

* 

r 

1 — r 

-1 


1 - 

1 _ 

J —00 



-1 

-1 


w{jw)_ 


(25) 


Thus Sr^r satisfies the IQC defined by the following multiplier: 


n4 := 


-1 

1—r 

-1 -1 


(26) 


This multiplier is analogous to the multiplier Hi given in 
Example 1 for constant delays. Hi corresponds to a unit 
circle in the Nyquist plane centered at —1. At each frequency 
n 4 can also be viewed as a circle centered at —1 but with 
radius enlarged to J— > 1 to account for the varying 
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delays. However this is not a precise statement since the 
constraints defined by IQC(n 4 ) only hold when integrated 
over all frequencies. 


Example 5. The multiplier in Example 4 depends on the rate 
of variation r but does not depend on the maximum delay f. 
Proposition 2 in fl^ gives a delay-dependent IQC that can be 
used to reduce the conservatism. The IQC in lfT9ll depends on 
a rational bounded transfer function (j) 5 {s) that satisfies; 


ksO'w)! > 


f\Lo\ ifT|u;|<l + ^7^ 

l + 7ib + 


(27) 


If r < 1 then Sf,r satisfies the IQC defined by the following 
multiplier Hs: 


n5(jw) 


l^5(jw)P 0 

0 -1 


(28) 


The proof that S IQC(n 5 ) uses the bound ||5r .r-O ill < 
r. The IQC multiplier IIs is analogous to the QC multiplier 
112 provided in Example 2 for constant delays. Specifically, 
a Taylor series expansion for Sf can be used to show that 
IIs is equivalent to 02 at low frequencies. The bound on 
effectively increases the radius of circle constraints defined 
by Hs at high frequencies to account for the time variations 
and to cover all delays in [0,7=]. Note that as r —>■ 0, the 
high frequency bound in Hs becomes \(j> 5 {juj)\ > 2. Thus as 
r —0, Hs converges to the multipliers 112 used to cover all 
constant delays in [0, f]. Again, these interpretations of Hs are 
imprecise and only meant to provide an intuitive interpretation. 
Proposition 3 in fl9l gives a similar IQC multiplier that is valid 
for r < 2. 


Example 6. Example 3 in Section IV-B provided a QC 
multiplier Ha for constant delays corresponding to a circle 
centered at The benefit of this multiplier is that it 

defined a smaller circle than the multipliers given in Examples 
1 and 2. This frequency domain intuition can be used to derive 
a new, related IQC for time-varying delays. The new IQC 
depends on a rational transfer function (j)Q{s) that satisfies; 


l^6(jw)| > 



if |r|w| <1-1- 
if |t|w| >1 -1- 



(29) 


It is shown in Appendix [A| that if r < 1 then Sf,r satisfies the 
IQC defined by the following multiplier Hg; 


n 6 (ja;) 


\M^)? - il'5r(jw)P \Sr(juj) 
\Sf{ju:) -1 


(30) 


The proof that Sf r & IQC(ng) uses the bound || (5^ r—^Sf)o 
i|| < if. Again, a Taylor series expansion for Sf can be used 
to show that Hg is equivalent to the analogous multiplier for 
constant delays Hg at low frequencies. The bound on |(^g| 
effectively increases the radius of circle constraints defined by 
Hg at high frequencies to account for the time variations and 
to cover all delays in [0,7=]. 


Each of the multipliers can be generalized for the case where 
Sf^r is MIMO using the idea of scalings already introduced 


in Section IV-A Sf^r is not a time-invariant system and hence 
frequency-dependent scalings cannot be used. However, the 
linearity of Sf^r can be used to show that constant matrix 
scalings can be introduced into the multiplier. Eor example, 
n 4 remains a valid IQC multiplier for Sf^r if it is generalized 
to include any matrix X > 0; 


n4 ;= 


-X 

1 — r 

-X -X 


(31) 


A formal proof that Sf^r S IQC(n 4 ) is given in Proposition 
1 of HD. Reference Cl derives additional IQCs for time- 
varying delays. In particular, frequency-dependent scalings can 
be introduced into the multipliers for varying delays but, in 
this case, a swapping lemma must be used to account for 
the the time-variations in the delay. This section does not 
intend to provide an exhaustive review of IQCs for time- 
varying delays. Instead the main purpose is to demonstrate 
the benefit of the frequency-domain intuition provided by QCs 
for constant delays. Example 6 is a new IQC derived for time- 
varying delays using this intuition and it should be possible 
to derive additional useful IQCs using this approach. 


V. Time Domain Stability Analysis 


The previous section defined frequency domain constraints 
that describe the input/output behavior of a time delay. These 
were specified as QCs for constant delays and IQCs for 
time-varying delays. This section shows that, under some 
mild technical conditions, these constraints have an equivalent 
time domain representation (Section |V-A[ ). The alternative 
time domain representation for QCs and IQCs is used to 
derive dissipation inequality based stability conditions for de¬ 


layed nonlinear and parameter varying systems (Sections V-B 
and[V^. 


A. Time Domain IQCs 

Let n be an QC or IQC multiplier that is a ratio¬ 
nal and uniformly bounded function of joj, i.e. H G 
]jL^i+'" 2 )x(mi+m 2 )^ previously QCs hold point- 

wise in frequency and hence they are also valid when inte¬ 
grated over frequency. In other words, if H is a valid QC 
multiplier for a system A then it is also a valid IQC multiplier 
for the same system. Thus it is sufficient to provide a time 
domain interpretation for IQCs. 

The time domain interpretation is based on factorizing the 
multiplier as H = where M = G 

'I' G ^ (mi-i-ms) restriction to rational, bounded mul¬ 
tipliers n ensures that such factorizations can be numerically 
computed via transfer function or state-space methods ll43l . 
1^ . f23i . Such factorizations are not unique and two specific 
factorizations are provided in Appendix using state-space 
methods. Eor a general factorization, T* is assumed to be stable 
but may be non-square (n^ ^ m) and possibly non-minimum 
phase. 





















Next, let (v, w) be a pair of signals that satisfy the IQC in 


( |23| l and define z{juj) := 
be written as; 


w(juj) 


z{jui)*Mz{joj)doj > 0 


/ z{t)'^Mz{t) dt > Q 

Jo 

where 2 ; is the output of the LTI system 'h: 

i>{t) = A^ip{t) + B^iv{t) + B^2w{t), V'(O) = 0 

z{t) = + D^iv{t) + D^ 2 w(t) 


(33) 


(34) 


Thus signals v € 00 ) and w G L™^[0,oo) satisfy 

the IQC defined by 11 = if and only if the filtered 

signal z = 4/ [ ^ ] satisfies the time domain constraint in 
(0. Similarly, a bounded, causal system A satisfies the 
IQC defined by 11 = if and only if ( [3^ holds for 

all V G [0, 00 ) and w = A(v). To simplify notation, 
A G IQC(n) will also be denoted by A S IQC('I',M). 
Fig. [^provides a graphical interpretation for this time-domain 
form of A G IQC(T',M). The input and output signals of 
A are filtered through T* and the output z satisfies the time- 


. Then the IQC can 


(32) 


By Parseval’s theorem 11451 . this frequency-domain inequality 
can be equivalently expressed in the time-domain as: 


It is important to note that the time domain constraint (331 
holds, in general, only over infinite time intervals. The term 
hard IQC was introduced in ll2^ referring to the following 
more restrictive property; A satisfies the IQC defined by 11 and 
Jq z{t)'^Mz{t) dt > 0 holds for all T > 0, z; G 00 ) 

and w = A(z;). By contrast, IQCs for which the time domain 
constraint need not hold over all finite time intervals are called 
soft IQCs. Hard and soft IQCs were later generalized in 


to include the effect of initial conditions and the terms were 
renamed complete and conditional IQCs, respectively. The 
hard/soft terminology will be used here. 

This distinction between hard and soft IQCs is important 
because the dissipation inequality theorems developed below 
require the use of time-domain constraints that hold over all 
finite-time intervals. One issue is that the factorization of 
n as T''^M'I' is not unique and hence there is ambiguity 
surrounding the hard/soft IQC terminology. As a result, the 
characterizations of hard and soft are not inherent to the IQC 
multiplier H but instead depend on the factorization (4/, M). A 
more precise definition that refers to the factorization (T', M) 
is now given. 

Definition 4. Let H be factorized as with T* stable. 

Then (T', M) is a hard IQC factorization of H if for any 
bounded, causal operator A G IQC(n) the following time- 
domain inequality holds 

rT 


domain inequality in (33i. A simple example is provided 


to illustrate the connection between the time domain and 
frequency domain constraints. 


/ 


z{t)^ Mz{t) dt > 0 


(37) 


for all T > 0, z; G L^^O.oo), w = A(w), and z = Tz ["]. 



Fig. 6. Graphical interpretation of the IQC defined by 11 = 


Example 1. Consider the multiplier n 2 (ja;) = 


I02(itj)| 

0 


-1 


(Equation 201 where (j )2 is a stable LTI system that satisfies 
certain bounds. This multiplier is a rational approximation 
that corresponds, at low frequencies, to a circle centered at 
the origin with a frequency-dependent radius. G QC(n 2 ) 
implies Sr G IQC(ri 2 ). Hence any v = L2[0, oo) and 
w = SrV satisfy 


v{jui) 

w(ju>) 


112 (jw) 


v(juj) 

w{ju>) 


dw > 0 


(35) 


A factorization for this multiplier is given by M = [ o _?i ] 
and 4! = «]. If zu = SrV then z := 4- [" ] = 

Hence the time-domain form for the IQC is 

pOO pOO 

/ z{t)'^Mz{t) dt — / zf{t) — zl{t)dt>0 (36) 

Jo Jo 


where Z 2 = w and zi = (l> 2 V. In other words, w has L 2 norm 
less than the filtered signal (j) 2 V, i.e. ||zz;|| < ||^i 2 t^||- 


It was recently shown that a broad class of IQC multipliers 
have a hard factorization ll24l . The proof uses a new min/max 
theorem to obtain a lower bound on z(t)'^Mz(t) dt. A 
similar factorization result was obtained in ll40ll using a game- 
theoretic interpretation. The next theorem summarizes the 
main factorization result needed in order to use IQCs within 
the dissipation inequality framework. 


Theorem 1. Let H = H'^ G 


titioned as 


•mi+m2)y.(mi+m2) 


^21 

1122 


par- 
and H,, G 


where Hu G iK.jLigQ* * aiiu 1122 

Assume nii(ja;) > 0 and n 22 (jw) < 0 for all 
a; G K U { 00 }. Then H has a hard factorization (T*, M). 


Proof. The sign definite conditions on Hu and n 22 ensure 
that n has a factorization (M, T*) where Tz is square and 
both are stable. This follows from Lemmas and 

in Appendix [B| Moreover, Appendix [B] provides a numerical 
algorithm to compute this special (J-spectral) factorization 
using state-space methods. The conclusion that (M, T') is a 
hard factorization follows from Theorem 2.4 in ll24l . □ 


B. Analysis of Nonlinear Delayed Systems 

This section derives analysis conditions for the nonlinear 
delayed system Fu{G,Sr) shown in Fig. [fusing dissipation 
inequalities. For concreteness the discussion focuses on Sr 
defined with constant delays but the results also hold using 
IQCs valid for Sf^r defined with varying delays. Assume Sr 
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satisfies the IQC defined by If and, in addition. If has a hard 
factorization M). The feedback system can be analyzed 
using the interconnection structure shown in Fig. with the 
system T* appended to the input/output channels of Sr- The 
dynamics of the interconnection in Fig. [^involve an extended 
system of the form 


X := F(x, w, d) 

= H (x, w, a) 

X := S R”< 3 +n,/, jjjg extended state and the functions 
F and H can be easily determined from the dynamics of 
G and T'. The theorem below provides a sufficient condition 
for the feedback interconnection to have an induced L 2 gain 
from d to e that is less than or equal to 7 . The theorem is 
based on a dissipation inequality that uses both the hard IQC 
associated with Sr as well as a storage function V defined on 
the extended state x. The system Sr is shown as a dashed box 
in Fig. [ 7 ] because the analysis essentially replaces the precise 
relation w = Sr{v) with the hard IQC constraint on z that 
specifies the signals pairs {v, w) that are consistent with the 
behavior of Sr- 


z 

e 



e 


Fig. 7. Analysis Interconnection Structure 


Theorem 2. Assume Fu{G,Sr) is well-posed and Sr satisfies 
the hard IQC defined by (T',M). Then ||F'„(G, 5.r)|| < 7 if 
there exists a scalar A > 0 and a continuously differentiable 
storage function V : 

i) 1 /( 0 ) = 0 , 

ii) vlx) > 0 Va; G 

iii) The following dissipation inequality holds for all x G 

Xz^Mz -G VI/(x) • F{x, w, d) < j^d^d — e^e (39) 

where z and e are functions of {x.,w,d) as defined by 
id in Equation [3^ 


Proof. Let d G L 2 ‘^[ 0 ,c») be any input signal. From 
well-posedness of the interconnection, the interconnection 
Fu{G,Sr) has a solution that satisfies the dynamics in Equa¬ 
tion The dissipation inequality (Equation [39| can be 


integrated from t = 0 to t = T with the initial condition 
x( 0 ) = 0 to yield; 


A f z{t)'^Mz{t) dt + V {x{T)) j d{tY'd{f)dt 
Jo Jo 


(40) 


— / e{t)^e{t)dt 

Jo 

It follows from the hard IQC condition, A > 0, and the non¬ 
negativity of the storage function V that 

f e{t)'^e{t) dt < f d{t)^d{t)dt (41) 
Jo Jo 

Hence |li^„(G,5,)|i < 7 . □ 


The dissipation inequality (Equation [39l l is an algebraic 
constraint on variables {x,w,d). The dissipation inequality 
only depends on Sr via the constraint on z'^Mz. Thus the 
dependence of the dissipation inequality on the delay value t 
appears through the choice of the multiplier II. Specifically, H 
typically depends on the value of t, e.g. II 2 and Tfa defined 
previously. The delay value is selected and then the multiplier 
n and its hard factorization (T*, M) are constructed. Thus for 
a given delay value t. Theorem provides convex conditions 
on V, A, and 7 that are sufficient to upper bound the Lo gain 

This leads to a useful numerical procedure if additional 
assumptions are made on G. If the dynamics of G (Equation]^ 
are described by polynomial vector fields then the functions 
F and F[ in the extended system are also polynomials. If the 
storage function V is also restricted to be polynomial then 
the dissipation inequality (Equation [39| and non-negativity 
condition 1 / > 0 are simply global polynomial constraints. 
In this case the search for a feasible storage function V and 
scalars A, 7 can be formulated as a sum-of-squares (SOS) 
optimization ||33|, ll22l . Eor fixed delay t this yields a 
convex optimization to upper bound the L2 gain of Fu{G, Sr)- 
In addition, bisection can be used to find the largest delay 
f such that the gain from d to e remains finite. If the QC 
multiplier H covers Sr for all r G [0, f] then f is a lower 
bound on the true delay margin. It is a lower bound because 
the dissipation inequality is only a sufficient condition. An 


example of this SOS method is given in Section VTA One 


issue with the SOS approach is that the required computation 
grows rapidly with the degree and number of variables con¬ 
tained in the polynomial constraint. This currently limits the 
proposed approach to situations where the extended system 
roughly involves a cubic vector field and state dimension less 
than 7 — 10. 

It should also be noted that the multiple IQCs can be 
incorporated in the analysis. Specifically, assume Sr satisfies 
the hard IQCs defined by for k = I,-- - ,N. 

Each T'fc can be appended to the inputs/outputs of Sr to 
yield a filtered output Zk- Theorem remains valid if the 
dissipation inequality (Equation [39) 1 is modified to include the 
term J2k=i ^kZ^MkZk for any constants Afc > 0. In this case 
the extended system includes the dynamics of G as well as the 
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dynamics of each 'I'fe (fc = 1, • • • ,N). The stability analysis 
consists of a search for the storage V, gain bound 7, and the 
constants Xk that lead to feasibility of the three conditions in 
Theorem This approach enables many IQCs for Sr to be 
incorporated into the analysis. 


C. Analysis of Parameter-Varying Delayed Systems 

A similar analysis condition can be derived for parameter- 
varying delayed systems. In particular, Linear Parameter Vary¬ 
ing (LPV) systems are a class of linear systems whose state 
space matrices depend on a time-varying parameter vector 
p : IR+ —. The parameter is assumed to be a continuously 
differentiable function of time and admissible trajectories 
are restricted, based on physical considerations, to a known 
compact subset V C K”'’. The state-space matrices of an 
LPV system are continuous functions of the parameter, e.g. 
Ac : V —> Define the LPV system Gp with inputs 

(w, d) and outputs {v, e) as: 


icii) = AG{p{t))xG{t) + Bcipit)) 


w(t) 

d(t) 


v(t) 

e(t) 


= CG{p{t))xG{t) + Dcipit)) 


w(t) 

d(t) 


(42) 


The state matrices at time t depend on the parameter vector 
at time t. Hence, LPV systems represent a special class of 
time-varying systems. Throughout this section the explicit 
dependence on t is occasionally suppressed to shorten the 
notation. 

By loop-shifting, a delayed LPV system can be modeled 
as Fu{G p,Sr) where w = Sr{v). This is similar to the 
interconnection shown in Figure but with Gp as the “nomi¬ 
nal” system. As a slight abuse of notation, ||F„(G'p, 5 t)|| will 
denote the worst-case L 2 gain over all allowable parameter 
trajectories: 

|jF„(Gp,5r)|| = sup sup (43) 

pGP 0^d<^L"'^[0,oo), xg(0)=0 


As in the previous subsection, assume Sr satisfies the IQC 
defined by H and that H has a hard factorization (fir, M). The 
stability and performance of Fu{Gp,Sr) can be assessed by 
appending the system T' to the input/output channels of Sr 
(as in Fig. |^. The dynamics of the interconnection in Fig. 
depend on an extended LPV system of the form: 


X = A{p)x Bi{p)w -f B 2 {p)d 
z = Gi{p)x + Dii{p)w + Di 2 {p)d (44) 

e = G 2 {p)x -V D 2 i{p)w -V D 22 {p)d, 


where the state vector is a; := ] S ^ 

denoting the state vectors of the LPV system Gp (Equation |42]i 
and the filter T' (Equation [34| i, respectively. A dissipation 
inequality can be formulated to upper bound the worst-case L 2 
gain of Fu{Gp, Sr) using the system (44 1 and the time domain 
IQC ( [J7| . This dissipation inequality is concretely expressed 
as a linear matrix inequality in the following theorem. The 
theorem is stated assuming a single multiplier H for Sr but 
many IQC multipliers can be included as described in the 
previous section. 


Theorem 3. Assume Fu{Gp, Sr) is well posed and Sr satisfies 
the hard IQC defined by (T'jM). Then ||Ftj(Gp, iSr)|| < 7 if 
there exists a scalar A > 0 and a matrix P = P'^ e 
such that P >0 and for all p G V: 


'A^P + PA 

PB, 

PB 2 ' 


'cf 

BfP 

0 

0 

+ 


1- 

to 

0 

1 

to 


.0^2. 


D 21 D 22 


+X 


Cl 

Dl, 

DI 2 . 


M 


Gi Dn D 


12 


< 0 


(45) 


In (451 the dependency of the state space matrices on p has 
been omitted to shorten the notation. 


Proof. Define a storage function V : 


IR+ by 

V{x) = x^Px. Left and right multiply (451 by [x’^, w'^, d'^] 
and to show that V satisfies the dissipation 

inequality: 


Xz{t)^Mz{t) + V{t) < 'Y‘^d{t)'^d{t) — e{t)'^e{f) (46) 


The remainder of the proof is similar to that given for 
Theorem |2l □ 


The analysis of the delayed LPV system reduces to a set of 
parameter dependent LMIs. These are infinite dimensional and 
hence they are typically approximated by finite-dimensional 
LMIs evaluated on a grid of parameter values. In this case the 
search for the matrix P and scalars A, 7 can be performed as a 
semidefinite programming optimization ||4l. If the LPV system 
has a rational dependence on the parameters then a finite 
dimensional LMI condition can be derived (with no approx¬ 
imation) using the techniques in ESll . ifTl . In addition, there 
may be known bounds on the parameter rates of variation. 
Theorem does not incorporate such knowledge and hence 
this is called a rate-unbounded analysis condition. Theorem 
can be easily extended to include rate-bounds using parameter- 
dependent storage functions as described in ll42l . Einally, if the 
system Gp is LTI then Equation 46 reduces to a single LMI 
condition. The nonlinear dissipation inequality in Theorem]^ 
is equivalent to exactly the same LMI condition when G is LTI 
and V is a quadratic function of x. In other words. Theorems]^ 
and are the same for LTI dynamics and quadratic storage 
functions. 


VI. Numerical Examples 

This section presents numerical examples to assess the 
stability and performance for nonlinear and LPV delayed 
systems. 




















11 


A. Nonlinear Delayed System 

Consider the classical feedback loop shown in Fig. where 
Dr is a constant delay. L is a nonlinear system described by 
the following ODE: 


XG = 

y = 


-49 0 

XG 

1 0 
-4.5 1.5 XG 


w+p{xg) 


(47) 


frequency and delay margin estimate of r = 1.96 sec (small 
red circle). The location of Sr lies inside the circle described 
by flop* at the critical frequency. In fact, Sr lies inside the 
circle described by li-opt at all frequencies. This is consistent 
with the following rough frequency-domain interpretation of 
the dissipation inequality result: Fu{G,A) is stable for LTI 
perturbations A whose Nyquist plot lies within the circles 
described Ilopt. 


where p{xg) := 2x‘^q+ Sx^ .^ - 0.2x1; -x^ ^ ■ 

described in Section]^ a loop-shift can be performed to bring 
the classical feedback loop into the form F„(G, Sr) shown in 

Fig-in 



Fig. 8. Classical Feedback Loop 


First, a linear analysis is performed to gain insight into the 
use of QCs for analysis. Let Gun denote the linearization 
of G around xg = 0. This linearization is obtained by 
neglecting the higher order terms in (47 1 , i.e. neglecting p(xg). 
An estimate for the delay margin of Fu{Giin,Sr) can be 
computed from the LMI condition in Theorem Bisection 
is used to hnd the largest time delay for which the gain 
from d to e is hnite. A number of solvers exist for this 
class of systems and here the Matlab LMI Lab toolbox was 
used. Using the standard multipliers IIi and 112 gives a delay 
margin of 0.06sec. However, using the multipliers Hi and the 
new “small” circle multiplier 143 yields a signihcantly larger 
delay margin estimate of 1.96sec. The exact delay margin for 
Fu{GiimSr) can be estimated from the frequency response 
of Gun- The frequency response gives a delay margin of 2.05 
sec with a critical frequency 0.361 rad/sec. These results are 
summarized in Tab. |I] 


System 

Method 

Delay Margin 

Linear, Gun 

IQC with Hi and n 2 

0.06sec 

Linear, Gun 

IQC with Hi and fts 

1.96sec 

Linear, Gun 

Frequency Response 

2.05 sec 

Nonlinear, G 

IQC with Hi and fta 

1.09sec 


TABLE I 

Summary of Delay Margins 


The results can be interpreted in the complex plane as 
shown in Fig. The hgure shows the circle constraints 
described by Hi and 113 the critical frequency 0.361 rad/sec 
(dotted blue). The hgure also shows the circle described by 
the optimal multiplier Ilopt '■= AiHi -f A 3 n 3 (solid blue). 
The coefficients in the optimal conic combination are given 
by Ai = 5310 and A 3 = 7210. Finally, the hgure shows 
the location of Sr{joj) = — 1 evaluated at the critical 



Real 


Fig. 9. Nyquist plot showing quadratic constraints. 


This frequency domain interpretation can also be used to 
explain the poor delay margin bound obtained with Hi and 
n 2 . Briehy, the transfer function from w to v (upper left 
comer of Gun) is given by —T = i.e. the negative 

complementary sensitivity function. Standard robust control 
techniques BTIl can be used to show that A := T(j^o) 
is a complex perturbation that causes Fu{Giin, A) to have 
poles at ijwo- At high frequencies these “small” destabilizing 
perturbations follow a vertical line in the complex plane with 
real part equal to 0.37. This constrains the use of 142 because 
this multiplier dehnes a circle centered at the origin with radius 
greater than Sr ■ TTie destabilizing perturbations lie within this 
circle unless the time delay is chosen to be small. This issue 
can be eliminated by frequency weighting the multiplier 142 , 


e.g. weighting each entry by 


s+r 


. However, this increases 


lOs-Fl 

the state order of the multiplier and hence the computation 
time. 

Next consider the analysis of the delayed nonlinear sys¬ 
tem. An estimate for the delay margin of Fu{G,Sr) can be 
computed from the dissipation inequality in Theorem A 
number of solvers exist for the corresponding sum-of-squares 
optimizations. Here the SOSOPT toolbox was used for all 
computations. Using the multipliers Hi and 143 yields a delay 
margin of = 1.09 sec for the nonlinear system. Thus the 
nonlinearities signihcantly degrade the delay margin bound. 
The bisection was initialized with the bounds tnl € [0,5]. 
The computation took 5.2 secs to perform 13 bisection steps 
to achieve a tolerance of 0 . 001 . 
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Fig. 10. Induced L 2 gain versus delay 


The L 2 gain of the delayed nonlinear system from reference 
d to tracking error e can also be computed from the dissipation 
inequality in Theorem Figure 10 shows the gain of the 
nonlinear system for delays below the delay margin estimate, 
T < tml = 109 sec (green dash-dot). The multipliers IIi and 
Tla were used to compute this curve. It took 7.5 sec to evaluate 
the gain on a grid of 20 delay values. For comparison the 
figure also shows the gain of the linear system Fu{Giin,Sr) 
computed using two methods. The red dashed curve is the gain 
computed using the LMI condition in Theorem also with 
multipliers IIi and Fla. The blue solid curve is the true induced 
L 2 gain of the linear system estimated from the frequency re¬ 
sponse of Fu{Giin, Sr)- The two linear results are close which 
provides confidence in the upper bounds computed via the IQC 
and dissipation inequality condition. The hgure also shows that 
the nonlinear delayed system has significantly larger gain as 
compared to the linearized system. This again indicates that 
the nonlinearities significantly degrade the performance. 


blade. Since Gp only depends affinely on the parameter p, it is 
sufficient to only consider the vertices, i.e. p = — 1 and p = 1 , 
for the following analyses. The goal of the benchmark is to hnd 
the time delay margin of the plant for a given cutting stiffness 
k. As in the previous example a loop shift is performed to 
bring the system into the form Fu{Gp,Sr)- 

The results of the analysis are shown in Fig. [TT] The IQC 
approach is based on the LMI condition Theoreml^and solved 
via bisection using the Matlab LMI Lab toolbox. The IQC 
analysis is performed using the combination of multipliers IIi 
and ri 2 - The IQC analysis is also performed using IIi and II 3 . 
Finally the method from ll44ll based on Lyapunov-Krasovskii 
functionals is shown for comparison with the IQC approach. 
It is shown in ll44l that the system is stable independent of the 
delay for stiffness values k < 0.267. As can be seen in the 
hgure, all three methods capture this behavior well. For higher 
stiffness values, the IQC method with using multipliers IIi and 
n 3 yields improved margin bounds as compared to the analysis 
condition from ll44ll . It should be noted that the improved 
Lyapunov-Krasovskii condition in ifTSll achieves results similar 
to the best IQC results obtained with IIi and II 3 . 



B. LPV Delayed System 

The second example is an LPV time-delayed system repre¬ 
senting a cutter used in milling. The system Gp, taken from 
ED, can be written as 


Fig. 11. Time Delay Margin vs Cutting Stiffness 


VII. Conclusions 


XG 


0 0 1 

0 0 0 

-(10-f 0.171fc)-f O.S/cp 10 0 

5 -15 0 


0 

0 


0.171fc- O.S/fcp 


0 


0 0 
0 0 
0 0 
0 0 


0 

0 

0 

0 


w 


W ='Dr{XG) 


0 

1 


-0.25 


(48) 


where k is the cutting stiffness and p € [— 1 , 1 ] is an artihcial 
scheduling parameter depending on the angular position of the 


This paper developed stability and performance analysis 
conditions for nonlinear and LPV time-delayed systems. The 
approach bounds the behavior of the time delay using integral 
quadratic constraints (IQCs). IQCs are typically specihed as 
frequency domain constraints on the inputs and outputs of an 
operator. For constant time delays, such IQCs were shown to 
have an intuitive geometric interpretation at each frequency. 
This intuition led to the construction of new IQCs for both 
constant and varying delays. Dissipation inequalities were then 
provided that incorporate IQCs into the analysis of delayed 
systems. This analysis approach applied existing results to 
obtain an equivalent time-domain intepretation for IQCs. Fi¬ 
nally, numerical examples were provided to demonstrate the 
proposed methods. While this paper only considers a single 
time delay, the framework can easily incorporate additional 
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delays and uncertainties analysis. Future work will explore the 
connections between the IQCs/dissipation inequality approach 
and standard Lyapunov methods. 
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Appendix 

A. Time-varying IQC 

This appendix provides a proof that if r < 1 then Sf,r € 
IQC(n6) where the multiplier Hg is defined in Section 
Two preliminary lemmas are required to prove this result. 

Lemma 2. If r < 1 then \\Sf^r — < 1 + yfry- 

Proof. Apply the triangle inequality as well as the definitions 
of Sf,r and Sf to obtain the following bound; 

- 2‘^rll - \\'^f,r\\ + 211^^11 + 2 


IV-C 


(49) 
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The varying delay is bounded as ||T’T,r|| < ifT^ . ifT^ 

while the constant delay is bounded by |jl?r|l < 1 - D 

Lemma 3. ||(5T,r — o ^|i < 

Proof. The proof is only sketched as it is similar to that given 
for Lemma 1 in ifTSl . To simplify notation define A := {Sr,r — 
iiSr)oi. Consider w = Au for some v G L 2 [ 0 , oo) and define 
y(t) := fg v(a)da. Thus w = {Sf^r — ^Sf){y) which, after 
some algebra, gives 

w{t) = f s{a)v{a)da (50) 

J t — T 

where s{a) = +\ for a G [t — f,t — T{t)\ and s(a) = — | for 
a G[t — T{f),t\. The Cauchy-Schwartz inequality can then be 
used to show 

— pt 

^ T / v‘^{a)da (51) 

4 Jt-i 

Integrate this inequality from f = 0 to f = c» and perform a 
change of variables to obtain Ikf < A||^|| 2 . □ 


Theorem 4. Let 4)^{s) be the transfer function defined in 
Equation 29 If r < 1 then |j(iSr,r — o < 1. 


Proof. The proof is only sketched as it is essentially the same 
as that given for Proposition 2 in 1191 . Let u G L 2 be an input 
signal and v := F{v) its corresponding Fourier Transform. 
Decompose u as ul + vh where vl and vh are the low 
and high frequency components, respectively. Specifically, the 
low-frequency content is defined in the frequency domain by 
VLijuj) ■■= v{juj) if |a;| < f + 7 ^) and VL{juj) := 0 
otherwise. The high-frequency content is defined similarly. To 
simplify notation, define A = {Sf,r — Then using 

the linearity of A and the triangle inequality yields 


||Au|| < IjAuill -f ||A'u/i-|| (52) 

Lemmas and bound the gains on the high and low fre¬ 
quency components by ||Ar!/i-|| < ||t;_fr|| and ||Aul|| < ||ul||. 
Thus ||Au|| < IIulII -f ||uh|| = ||u||. □ 


The bound in Theorem can be equivalently expressed as 
a quadratic, frequency-domain constraint on the input/output 
signals of Sf^r- This gives the desired result that Sf r satisfies 
the IQC defined by the multiplier fig. 


B. IQC Factorizations 

This appendix provides specific numerical procedures to 
factorize 11 = 11'" G as 'I''"MT'. Such factoriza¬ 

tions are not unique and this appendix presents two useful 
factorizations. 

First, let be a minimal state-space real¬ 

ization for n. Separate 11 into its stable and unstable parts 
n = Gs + Gu- Let ( 2 I, B, G, i9^) denote a state space 
realization for the stable part Gs. The matrix A is Hurwitz 
since Gs is stable. The assumptions on 11 can be used to show 
that the poles of 11 are symmetric about the imaginary axis 


and, moreover, Gu has a state space realization of the form 
(-A^, -G'^, B'^, 0) (Section 7.3 of lIH). Thus H = Gs+Gu 
can be written in the form 11 = where 


T'(s) := 


M := 


{sI-A)-^B 

I 

0 G^' 

G 


(53) 


(54) 


This provides a factorization 11 = where M = M'^ G 

and T* G poj- this factorization T* is, in 

general, non-square (n^ f m) and it may have right-half plane 
zeros. 

The stability theorems in this paper require a special factor¬ 
ization such that T* is square {nz = m), stable, and minimum 
phase. More precisely, given non-negative integers p and q, let 

. T* is called a Jo 


Jp^q denote the signature matrix 


Ip 0 
0 


p,q 

mxm 

•00 


spectral factor of 11 if 11 = T''" Jp and T', G 
The term J-spectral factor will be used if the values of p 
and q are not important. J-spectral factorizations have been 
used to construct (sub-optimal) solutions to the Hao optimal 
control problem m, Eol, im. The next lemma provides a 
necessary and sufficient condition for constructing a J-spectral 
factorization of 11 . 


Lemma 4. Let 11 G 


n(s) = 


(sI-A)-^B 
I 


be a multiplier in the form: 

{sI-A)-^B 
I 



0 

G^' 



G 




(55) 


where A is Hurwitz. The following statements are equivalent: 

1 ) Dtt = is nonsingular and there exists a unique real 
solution X = X'^ to the following ARE 

A'^X -f AA - {XB + G^)Df^{B^X -f C) = 0 

(56) 

such that A - BDf^ {B^X + G) is Hurwitz. 

2 ) n has a Jp ^ spectral factorization where p and q are 
the number of positive and negative eigenvalues of 
Dtt, respectively. Moreover, 'L is a Jp ^-spectral factor 
of n if and only if it has a state-space realization 
[A, B, Jp^qW~'^ [B"^X 4-G) ,W) where VF is a so¬ 
lution of = W'^ Jp^qW. 

Proof. This lemma is based on the canonical factorization 
results in lO and summarized in Chapter 7 of ifTTIl . The precise 
wording of this lemma is a special case of Theorem 2.4 in 

ca. □ 


The numerical procedure to construct a J-spectral factoriza¬ 
tion of n = H'" can be summarized by the following steps. 
First, express H with a minimal realization (A^, B^^, D^f). 

Second, compute a state space realization {A,B,C,Dpf} for 
the stable part of H. This step can be done using the Matlab 
command stabsep. Third, attempt to solve for the stabilizing 
solution X = of the ARE in Equation 56 This step can 
be done using the Matlab command care. The existence of a 
stabilizing solution is related to the eigenvalues/eigenspaces 
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of a related Hamiltonian matrix (see ||45]| '). If this step is 
unsuccessful, i.e. no such solution exists, then H does not have 
a J-spectral factorization by Lemma However, if the ARE 
has a unique stabilizing solution then construct the state-space 
realization for as defined in Statement 2) of Lemma This 
requires the matrix decomposition = W"’'Jp^qW which 
can be computed from an eigenvalue decomposition of 
This entire procedure, if successful, yields the factorization 
n = where M := and T*, S 

The last result of this appendix provides a simple frequency 
domain condition that is sufficient for the existence of a J- 
spectral factor of a multiplier H = H'^. 


Lemma 5. 
tioned as 


Let n = n 

till ni2 

n~2 n22 


e 

where Hu 


mi+m2)x(mi-rm2) 

: and H,, G 


Assume nii(ja;) > 0 and n 22 (jw) < 0 for all 
(jj G ffi U {oo}. Then H has a Jmi,m 2 ‘SP 6 ctral factorization. 


Proof. The sign dehnite conditions on Hu and n 22 can be 
used to show that H has no equalizing vectors (as dehned 
in Ezl) and hence the corresponding ARE has a unique 
stabilizing solution (Theorem 2.4 in 1271 1. Details are given 
in EOl. □ 




